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ABSTRACT 

We have investigated an influence of the r-mode instabihty on hypercriticaUy accreting 
(M ~ lM0y~^) neutron stars in close binary system s duri ng their common envelope 
phases based on the scenario proposed by Bethe et al. ( |l999| ). On the one hand neutron 
stars are heated by the accreted matter at the stellar surface, but on the other hand 
they are also cooled down by the neutrino radiation. At the same time, the accreted 
matter transports its angular momentum and mass to the star. We have studied the 
evolution of the stellar mass, temperature and rotational frequency. 

The gravitational-wave-driven instability of the r-mode oscillation strongly sup- 
presses spinning-up of the star, whose final rotational frequency is well below the 
mass-shedding limit, typically as small as 10% of that of the mass-shedding state. On 
a very short time scale the rotational frequency tends to approach a certain constant 
value and saturates there as far as the amount of the accreted mass does not exceed 
a certain limit to collapse to a black hole. This implies that the similar mechanism of 
gravitational radiation as the so-called Wagoner star may work in this process. The 
star is spun up by accretion until the angular momentum loss by gravitational radi- 
ation balances the accretion torque. The time-integrated dimensionless strain of the 
radiated gravitational wave may be large enough to be detectable by the gravitational 
wave detectors such as LIGO II. 
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1 INTRODUCTION 

In relation to the formation mechanism of compact binary 
systems such as double neutron star systems, double black 
hole systems or black hole-neutron star systems, hypercriti- 
cal accretion flows onto neutron stars in envelopes of massive 
stars of close binary systems have been invest i Rated by sev- 



eral authors (Chevalier 1993 



Chevalier 1996 



Brown 1995 



Bethe fc Brown 199^ ; pethe et al. 19991 ). The mass accre- 



tion rate {M) of the hypercritical accretion flow amounts to 
10* times as hi^ h as the Eddington's mass accretion limit. 
Chevalier (19961) found that for rather small values of the vis- 
cosity parameter, a ~ 10~^, the accretion flow which cools 
via neutrino radiation coul d be r ealized. 

Recently Bethe et al. ( |1999| ) pointed out that the flow 
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with much larger and conservative values of the viscosity pa- 
rameter, a ^ 0.1, does not become so hot as to emit neutri- 
nos during the accretion process but that the accretion will 
proceed to the advection dominated accretion flow (ADAF). 
In their model, the accreted matter onto the neutron star 
surface is processed there by nuclear reactions and cools via 
neutrino emission. The hypercritical accretion increases the 
mass of the star rapidly and induces the collapse of the neu- 
tron star to a black hole. Thus the supernova explosion of 
the massive companion will lead to the formation of a black 
hole-neutron star binary system or a double black hole bi- 
nary system. The merging event of such a system might be 
one of the important targets for the gravitational wave as- 
tronomy. 

In the accretion process, a non-negligible amount of an- 
gular momentum is carried into the neutron star. This po- 
tentially could spin up the star to the Kepler limit, i.e. to the 
state where the stellar rotational frequency equals to the or- 
bital frequency at the equatorial surface. The star rotating 
at the Kepler limit deforms significantly from the spheri- 
cal configuration and the subsequent collapse of the rapidly 
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rotating star to a black hole might radiate a substantial 
fraction of the binding energy of the system in the form of 
gravitational waves. 

This seemingly straightforward conclusion, however, 
has to be reexamined from the standpoint of stability of 
rotating neutron stars. The gravitational radiation driven 
instability of rotating stars might prevent the stellar rota- 
tional frequency to increase above a certain limit . In partic- 
ular , the recently discovered r-m ode instability (Andersson 
1995 ; F iedman & Morsink 1998) might be sufficiently strong 



to suppress the spin frequency in some range of the stellar 
temperature. 

In this paper, in order to see how this instability influ- 
ences the accretion induced spin-up of the neutron star in 
the hypercritical accretion, we will employ the hyperc ritica l 
ADAF model based on the analysis of Bethe et al. ( 199£ ) 
and investigate the evolution of the temperature and the 
rotational frequency of the neutron star. ^ 



2 SIMPLIFIED MODEL OF THE ACCRETION 
INDUCED EVOLUTION OF NEUTRON 
STARS 

It should be noted that the study in this paper is done 
in the framework of Newtonian dynamics and Newtonian 
gravity. The gravitational radiation reaction on stellar os- 
cillations is eva luated by the lowest order post-Newtonian 
approximation (Thorne 1980). Furthermore to simplify our 
analysis, stellar configurations are assumed to be spherical 
polytropes, p — Kp^, with the polytropic index A*' = 1. Here 
p and p are the pressure and the density, respectively. By 
choosing the polytropic coefficient K appropriately, the stel- 
lar mass is M = IAMq and its radius is i? = 12.5km. Note 
that the radius of the polytropes with index A*' = 1 is inde- 
pendent of the central density and thus of the mass if the 
constant K is fixed. 



2.1 Characteristics of the hypercritical ADAF 
onto neutron stars 

The hypercritical accretion flow whic h we employ in this pa- 
per is proposed by Bethe et al. (1999). The characteristics of 
that flow are summarized as follows. (1) The flow is hyper- 
critical because the accretion rate is typical ly ^ IM^y" ^ 



(2) S ince the flow is advection dominated (Narayan & Yi 



1994) and its temperature is relatively low (~ 0.5MeV), the 
adiabatic exponent of the gas equals to 4/3 (radiation dom- 
inated mixture of photon, electron, positron and nuclei). (3) 
The accreted matter touches down on the surface of the star 
'softly' and burns to produce as heavy nuclei as ^'^Ni. On the 
other hand, the matter is cooled by neutrino emission. Con- 
sequently the whole process acts as a 'thermostat' which 
keeps the neutrino temperature ~ IMeV. 



3 Similar investigations of accreting neutron stars with much 
lower accretion rate, which are the models of low mass X -rav 
binaries, can t^ e seen for instance in Andersson et al. ( 1999| ) and 
in Levin (|l99s|). 



2.2 Mass accretion rate 

The mass accretion rate A/ac 
attachment to the star, is assumed to be constant during 
the evolution. In general the mass accretion rate should be 
distinguished from the mass inflow rate Min, which is the 
total amount of the incoming mass per unit time to the 
neutron star. This is because some kinds of outflow may 
coexist with the accreting flows. In this paper we simplify 
the situation by assuming that both rates are the same, i.e. 
Afin = A/acc = M. We will discuss this assumption further 
in Discussion. 

The evolution of the neutron star mass M is written as: 



M 



M, 







^0 + /3t , 



(1) 



where is the initial mass of the star (in units of the solar 
mass) , time t is measured in units of second, and /3 is deflned 
as; 



/3 = 3.2 X 10" 



M 



IMey- 



(2) 



2.3 Thermal response of the neutron star heated 
by neutrino from the surface 

The accreted matter landing on the surface layer is processed 
by nuclear burning and emits neutrinos. We will assume that 
roughly a half of them escapes from the star freely and that 
some fraction of the rest of them interacts with stellar matter 
and heats it up. The mean free path A of inelastic scattering 
of neutrinos with electrons in the neutron star matter is 



written as (Shapiro & Teukolsky 1983), 



A ~ 2 X 10' 



Pnuc 


7/6 


rO.lMeVi 


5/2 

cm 


p 




[ \ 





(3) 



where pnuc is the 'nuclear density' (~ 2.8 x 10 g cm~ ) and 
is the neutrino energy. The density p means the averaged 
density of the star, which is related in our polytropic stellar 
models to the mass M as follows: 



P 

pnuc 



: 0.7 



M 



(4) 



The neutrino energy is fixed to IMeV. We simply as- 
sume the scattering (thus heating) rate as r/R/X where R is 
the stellar radius and rj is an efficiency factor of f-heating 
through the inelastic scattering. As the flow is adiabatic and 
radiation dominated, the specific internal energy of the ac- 
creted matter s is written as. 



3 6GM 

4 7~R' 



(5) 



With these expressions, the heating rate of the whole star 
e'h is. 



X 



4.7 X 10* 



M 



lAfoy- 



M 



-1 13/6 



erg s 



(6) 



Here we choose r) = 1. Heating rate above is proportional to 
the multiplication of M and r). Thus the reduction of rj is 
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equivalent to that of M with rj being fixed. As shown later 
in Discussion, the dependence of the results on -q is weak. |^ 
As the cooling mechanism of the st ar, we assume the 



modified URCA process to be dominant ( Shapiro & Teukol- 
sky iggph. Its cooling rate L]^^'^^ is. 



6 X 10'' 



M 

Mr, 



2/3 



(7) 



where Tg is the stellar temperature, which is assumed to 
be uniform, in units of lO^K. The heat capacity Cv of the 
neutron star which is approximated by the following formula 
if the neutron star is assumed to consist of degenerate free 
nucleons, 



: 8 X X 



M 



Me 



-2/3 



(8) 



where fcs is the Boltzmann constant. Then the equation of 
the thermal balance of the star can be written as. 



(9) 



This is transformed to the equation of the evolution of the 
temperature as follows: 



di 



' M ' 


1/3 \ 






Me_ 





= 6.5 X 10" 



.7 X 10" 



M 



IMqy 
M 



M 



Mp 



Mp 



TS 



13/6 



(10) 



2.4 Angular momentum balance of the star 

The angular momentum of the star is changed by two 
torques induced by accretion and by gravitational radiation 
if the r-mode oscillation is unstable. 

The accretion torque is estimated as Mja, where ja is 
the specific angular momentum of the accreted matter: 



Ja 



R 



2 GM 
7~ 



(11) 



Here the ADAF solution by Narayan & Yi (1994) is em- 
ployed to evaluate the rotational velocity of the accreted 
matter at the stellar surface. 

As for the gravitational radiation torque, we use the 
simple formula used in Lindblom et al. ( 1998| ) which also in- 
corporates the damping effect by viscosity. This corresponds 
to the assumption that the r-mode quickly grows to saturate 
in the non-linear regime of its amplitude and the angular mo- 
mentum of the mode is a substantial fraction of the stellar 
angular momentum. 

Let us introduce the rate 5 of the angular momentum 
of the mode to the stellar angular momentum. By using the 
moment of inertia of the star I and the stellar rotational 
angular frequency fl, the torque due to gravitational wave 
emission can be written as rr^IQS. The factor S is related 



^ Here heating of the star by the r-mode oscillation, which is 
considered in Levin (1999), is not taken into account. We con- 
firmed, however, that an inclusion of the r-mode heating term is 
negligible in the situation of our interest. 



to the dimensionless amplitude of the oscillation of the r- 
mode 'k', which is used in the simplified discussion of the 
non- linear evolution of the r-mode by Owen et al. (1998), as 

The time scale of the r-mode instability Ti is then de- 
fined as. 



-1 _ -1 . -1 , -1 

— Terr I '''s + '''h 



(12) 



where rgr is the energy dissipation time scale by the grav- 
itational radiation, Ts is that of the shear viscosity due to 
neutron-neutron collision and rb is that of the bulk viscosity 
due to the lag of the /3— reaction and the stellar oscillation. 
The negative value of signals the onset of the instability. 
We set to be zero when this value is positive, i.e. for 
damping oscillations. 

The time scale of gravitational wave dissipation is ob- 
tained as, 



327rG {I - if 



1 + 2 

c2i+3 [(2/-f 1)!!]2 Vr+T 
,"1 



21+2 



p{x)x dx. 



(13) 



Here p is the density distribution normalized by the central 
density pc- The central density is related to the stellar mass 
M for N = 1 polytropes as. 



Pc = 3.29p = 3.29 X pnuc X 0.7^ 



(14) 



Thus for I — m — 2 modes for which the instability is the 
strongest. 



r"^ = -3.1 X 10" 



8 


r / 1 


6 


■ M ' 




LlOOHzJ 




Mq 



(15) 



Here / is the rotational frequency of the star (/ = f2/27r). 

The damping time scales due to shear and bulk viscosi- 
ties are, 

= {l-l){2l + l) / rjr^'^dr/ / pr^'+^dr 
Jo Jo 

„9/4p5Tn-lB^-2 r p(^f/\*dx 



5 X 347p^'"i?°10-'*r9"' 



2.0 X 10" 







pix)x dx 

M ^ 

Mo 



To 



(16) 



and. 



= (6.99X 10«)-^ f^Vrl 
\ttGpJ 



= 1.4 X 10" 



11 


r / 1 


2 


■ M ' 




LlOOHzJ 




Mq 



(17) 



The equation for the evolution of the rotational fre- 
quency is written as. 



d_ 

di 



M 



Mq 

+ 2.0 X lO"** 



lOOHz 

M 



IMqj- 



M 



Mc- 



M 

m; 

1/2 



f 



lOOHz 



(18) 
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2.5 Terminal values of the temperature and the 
rotational frequency of nearly constant mass 
models 

It is important to note tliat tlie tiiermal balance as well as 
the angular momentum balance determines the approximate 
terminal values of the temperature and the rotational fre- 
quency of the star, if the accretion rate is fixed and the mass 
is regarded as nearly constant in the evolution. 

By setting the right hand side of Eq. (^|) to be zero, 
the temperature can be solved as. 



Tg = 3.0 



M 


1/8 


■ M ' 


IMoy-i 







-1 7/48 



(19) 



Thus for the star with mass M ~ IAMq and the accretion 
rate with M ~ lAf0y~^, the (uniform) temperature will 
be about 3.5 x lO^K, which is near the most s usceptible 
temperatur e for the r-mode instability to grow (Lindblom 
et al. 1998| ). 

The equation of the angular momentum balance gives 
the terminal value of the rotational frequency of the star 
by setting the right hand side of Eq.(^^ to be zero. If we 
omit the contributions of viscous damping whose effect is 
smallest near at the temperature obtained above, we get 
the expression for the terminal rotational frequency as ; |IIJ 



r f 1 


= 1.3 


M 


1/7 


■ M ' 


LlOOHzJ 


lMQy-\ 







-3/14 



-1/7 



(20) 




1.0E6 2.0E6 



1.55 
Mass / Wau^ 



Figure 1. Evolutionary curves of tlie rotational frequency 
//lOOHz for different values of the initial rotational frequency, 
i.e. //lOOHz = 0.1,0.5,1 and 3, are plotted against the stellar 
mass which can be regarded as the 'time' variable. Model pa- 
rameters are as follows: the initial mass M = IAMq, the initial 
temperature Tg = 10"'', the mass accretion rate M = lMQy~^. 
Numbers attached to the curves are the values of the lapse of time 
in units of second. Also plotted is the curve of the model in which 
the r-mode instability is not included (dashed curve). 



2.6 Evolutionary curves 



Given the mass accretion rate and the initial values of the 
stellar mass, temperature and frequency, the evolution of 
the corresponding quantities can be followed by integrating 
Eqs. (|l]), (|l^) and (p^. We have computed several evolu- 
tionary sequences by choosing different sets for the initial 
values. 

We will show how the evolution depends on the initial 
values. In Fig. |^ the rotational frequency is plotted against 
the mass of the star which serves as the 'time' variable be- 
cause the mass accretion rate is assumed to be constant. 
With the mass accretion rate fixed, solutions with different 
values of the initial frequency converges to a certain value 
rather quickly compared with the mass accretion time scale. 
This behaviour is expected from the rapid temperature sat- 
uration seen in Fig. H by swift establishment of the thermal 
balance between heating and cooling. It should be pointed 
out that the omission of the r-mode instability would lead to 
very rapid increase of the rotational frequency and that the 
star would reach its mass-shedding limit before the accretion 
induced collapse occurs. 

The mass accretion rate dependence of the frequency 
evolution is shown in Fig. ^. As seen in Fig. ^ the tempera- 
ture of the star quickly saturates at the limiting value which 
depends on the mass accretion rate. Since the strength of 



1.55 
Mass / Mjun 



Figure 2. Evolutionary curves of the temperature Tg = 
T/10^(K) for different values of the mass accretion rate, i.e. 
M/(lMoy-l) = 1 (solid curve), 10'^ (dotted curve), IQ-^ 
(dashed curve), are plotted against the stellar mass. Model pa- 

1.4M0, the ini- 



rameters are as follows: the initial mass M 



tial temperature Tg 
//lOOHz = 0.5. 



10- 



and the initial rotational frequency 



II After the submission, we fin d a naner on an effect of a magnetic 
field to the r-mode instability ( Rezzolla et al. 200C ). Presence of 
the magnetic field may weaken the instability significantly. If this 
is the case with the situation of our interest here, the terminal 
frequency will be much higher. 
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Figure 3. The rotational frequency is plotted against the stellar 
mass. Three lines correspond to the three temperature lines in 
Fig. 2 (M/(lM0y-i) = 1 (solid curve), 10'^ (dotted curve), 
10~^ (dashed curve)). Numbers attached to these lines are the 
values of the lapse of time in units of second. 



the r-mode instability depends on the temperature as well as 
the stellar rotational frequency, the rotational frequency also 
quickly saturates. As the time scale of the angular momen- 
tum change is longer than that of the temperature change, 
the frequency saturation is reached slightly later compared 
with the temperature saturation. 



rate is considerably different, i.e. increased. The standard 
formula of th e luminosity of gravitationa l radiation can be 
expressed as (Shapiro & Teukolsky 1983), 



^^^h^D' ~ 2n fMja , (21) 

where uj,h,D are the eigenfrequency of the mode, the di- 
mensionless strain of the gravitational field and the distance 
from the source, respectively. Using the terminal rotational 
frequency / obtained above and the relation of the eigenfre- 
quency and the angular velocity for I = m — 2 modes, 

a; = i ■ 27r/ (22) 
we get 



5.6x10" 



22 


M 


3/7 


■ M ' 


5/14 


D 








Mq 




Ikpc 



cl/14 



(23) 



As seen from the result obtained in this paper, the rota- 
tional frequency becomes almost constant after a short time 
of accretion. Thus it is rather easy to detect signals of the 
gravitational wave by integrating it for an appropriate in- 
terval of time. 

We here define the Fourier component of the squared 
dimensionless strain as h? . Then from the Parseval's re- 
lation. 



4% 

-i* cycle 5 



(24) 



where A^cycio is the number of cycles of the oscillation in 
the interval of observation, Atobs. Frequency width Alj and 
Atobs are related by the uncertainty principle as. 



3 DISCUSSIONS 

3.1 Formation of slowly rotating black holes 

We have shown that the r-mode instability sets a severe 
limit on the rotational frequency of the hypercri ticall y ac- 
creting neutron stars proposed by Bethe et al. (199!:). As 



the mass accretion rate dependence of the rotational fre- 
quency (Eq. (^o|)) is rather weak, the rotational frequency 
of the star is well below its mass-shedding limit. Even if the 
mass accretion rate rises to M = lO^M0y~^, the rotational 
frequency increases only to 200Hz. Thus the gravitational 
collapse of neutron stars due to mass accretion to low mass 
black holes in these systems will proceed almost spherically 
and black holes formed just after the collapse will rotate so 
slowly that the rotational parameter is not large either, i.e. 
q = cJ/GM^ ~ 0.01 for 1.7Mq. 

3.2 Observational possibility of gravitational 
waves from the hypercritically accreting 
neutron stars 

The amount of gravitational radiation from the system can 
be evaluated as follows. The process d iscuss ed in this paper 
is similar to that studied by Wagoner ( 1984 ) . A neutron star 



in an X-ray binary system emits gravitational wave so as to 
balance the accretion torque with the gravitational radiation 
torque of the f-mode ('Wagoner star'). For our models the 
f-mode must be replaced by the r-mode and the accretion 



AlO ■ Atohs > - . 



(25) 



As Atobs = -^Ncycic, the r.m.s. of the Fourier component of 
the strain is. 



- V Atobs , 
h < -y- — h. 

- 21/4 



(26) 



The event rate TZ of the hypercritical accretion is 
quite uncertain. The maximum value is expected to be 



10" 



y 



"^galaxy ^, which is t he birth rate of the c lose 



black hole-neutron star binary (Bethe & Brown 1998). In 
order to detect A'o events per year, we have to observe the 
region with the distance as follows: 



D = (3iVo/47r7^ngal)'''^ 



(27) 



where ngai ~ O.OlMpc"'^ (Efstathiou et al. 1988). Here rigai 
is the number density of the massive galaxies (~ 10^^ M©). 
Then Fourier component of the strain is. 



/i(Hz 



-1/2. 



^ 1 X 10" 



rl/14 



M 



M 



-1 5/14 



Mc- 



IMoy- 

1/3 



lO'^galaxy 



-1 3/7 



-1/3 



to 



(28) 



1/2 



1 month 



Comparing this with the expected sensitivity of the Laser 
Interferometer Gravitational- Wave Observatory (LIGO), 
the integration of one month is sufficient for the detection 
by the LIGO II facility (see Shoemaker 2000 for the sensi- 
tivity curves of the LIGO). As the dependence of the strain 
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on TZ is weak, the LIGO II can see the signal by the inte- 
gration of one month, even if TZ is reduced by two orders of 
magnitude. Typical frequency of the source may be around 
100 ~ 200Hz. 



3.3 Uncertainties in our analysis 

In this investigation we have applied the 'one zone' approx- 
imation for simplicity. It means that the variables such as 
the stellar temperature or the rotational frequency are rep- 
resented by some 'averaged' values throughout the whole 
star. It premises that the viscosity works very efflciently 
compared to the time scale of the evolution of the system, 
typically about one year. The adjustment time scale of shear 
viscosity to smooth out the differential flow can be estirnated 



for the neutron-n eutron scattering dominant matter ( Lind- 
blom et al. 1998 ) as follows: 



-Tadj 



.Mm 



i/4- 



(29) 



Assumi ng the sLellai ladiuh Til — 10k ill cLil d Lilt; a V waged deu- 
sity p = 3 X 10"g/cm^, this amounts to 0.8 y for T = 10**K. 



The rol ation law of the star can deviate from the uniform 
rotation it this time scale is comparable or longer than the 
spin-up time. In that case, behaviour of (classical) r-modes 
is uncertain and their characteristics may be drastically 
changed. We here assume that they can be approximated 
by some kind of average of the rotational angular velocity. 
The validity of this assumption should be confirmed else- 
where. 

It has been pointed out that in the hypercritical ac- 
creti on process, explos ive outflows driven by neutrino heat- 
ing (Fryer et al. 199^ or jet formation (Armitage & Livic 



199!; ) may significantly reduce the mass attachment onto the 



neutron star. If it is the case, it is more probable that dou- 
ble neutron star binary systems will be formed through this 
channel. In such accretion inefficient situations, the mass ac- 
cretion rate Mace is lower than the mass inflow rate Min from 
the envelope of the companion star. If the accretion ineffi- 
cient fiows behave similarly as the accretion efficient flows, 
evolution of the stars may be approximated by our models 
with lower M and with a shorter duration, for the same value 
of Min- This results in smaller terminal rotational frequen- 
cies of the star. Thus neutron stars which survive after mass 
accretion in this process would have even lower rotational 
frequencies.^ 

It is a surprise that the possibility of observation of 
gravitational waves from the systems seems unexpectedly 
large. This might encourage the ongoing detection experi- 
ments of gravitational waves. However we have to be careful 
in concluding the significance of the hypercritically accreting 
neutron stars as sources of gravitational waves. Our analy- 
sis here might be too simplified and naive to get conclu- 
sive answers to the observational possibility. Some simpli- 
fied assumptions may influence the result significantly. For 



instance, more realistic behaviour of the accretion flow and 
the accreted matter may drastically change the efficiency of 
the accretion torque. Angular momentum re-distribution in 
the realistic star may modify the r-mode characteristics sig- 
nificantly, which may affect the strength of the instability. 
Thus the uncertainty included in our assumptions must be 
investigated more intensively in the future. 
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** It is possible that the characteristics of the accretion flow 
change drastically and that the accreted matter carries less an- 
gular momentum than the case considered here. Of course this 
will lead to the even more reduction of the terminal rotational 
frequency. 
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